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MICROHATO SCAITERING PRORERTIES OF SNOW FIELDS 
D. J. Ai'igclakos 

Electronic'S n.escarch Laboratory', 
University of California 
Berkeley, CA 947 2.0 

ABSTR^.CT 


A review is presented of the results obtained in tke research progretm of 
loicrovave sc.attering properties of snow fields. The program is a cooperative effort 
between a group at the University of California, Berke.ley (D. J. Angelakos, F. D. 
Clapp, K. K. Hex, and S. Coen) and NASA-Ames (W. I. Linlor) . A goodly portion of 
the research was undertaken at the Central Sierra Siiov; Lab and with the active 
cooperation of the U.S. Forest Service (J. L. Smith). 

Er.pei*imental results are presented showing backscatter dependence on a) 
frequency (5.S-8.Q CHz) , b) angle of incidence (0-60 degrees), c) snow v/etness • 
(time of day) , and d) frequency modulation (0-500 I'lHz) . Theoretical studies are 
being made of the inverse scattering problem yielding some prelirainary results 
concerning the determination of the dielectric constant of the snow layer. 

The experimental results lead to the following conclusions: a) Snow layerinj^, 

affects backscatter, b) Layer response is significant up to 45 degrees of incidence, 
c) Wetness modifies snow layer effects, and d) Frequency modulation "maslcs'* the 
layer response. 


latroduction 


At the Electronics Research Laboratory, University of California, Berkeley, the 
research p7:ogram of micrev/ave scattering properties of snow fields along two 
directions. One is the experimental detennination of the microvrave properties of ^ 
snow fields and the other is a theoretical program of mathematical modeling and the | 

CO 

inverse scattering problem. The ultimate goal of the project is the coalescence of ^ 
the two into, a useful diagrtostic technique for the determination of water content 
in snow . 


Experimental Portion of the Project - D. J. Angelakos, F. D. Clapp, W. I. Linlor 
and J. L. Smith: ITxe primary objeeti.ve of the experimental program was to deteT-mine 

the effect of snow conditions on the amount of backscattered e-lectroioagnetic energy. 
In addition the investigation, was to seek out any effects on the backscatter due to 
snow layering and to determine the dependence on frequency and frequency deviation. 

Ini proceeding towards the objective several additional phenomenological effects 
appeared which added insight to backscatter signal response. What follows is a 


sampling of results v^liich indicates the partial attainment of 
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Equipment ; Most of the data were taken for sv?ept-f requency ranging from 5.8 to 
I 8.0 KIz. Tlie sweepers the "Alfred" and the "Hewlett-Paclcard" sources. These 

iwere connected v;ith coax cable to horns having nominally 20-dB gain. The receiver 
jhqriis were coupled to crystal detectors and HP SWR meters, v^hose output v?as plotted 
'on. XY plotters, on a logarithmic scctle. For some measurements a converter unit was 
employed (Wiltron manufacturer) that permitted a uniform dB scale for the Y-axis. 

I For the 13.5 GH.z measurements a fixed-frequency solid state oscillator was 
used (Varian Eianufacturer) . Time \^s not available to take measurements at other 

i 

; frequencies, either fijied or swept. 

j Sites ; Two sites were used: one is called the "tow^er site" \diich is on the roof of 

jthe snow hut, for which only two fixed angles (39“ and 45®) could be used; the other 
|±s: called the "ladder site" where all desired angles could be obtained. 

The limited available time — - dictated by the snow depth — for measurements 

i 

|and also the limited manpower made a choice necessary between taking data or ground 
’truth (i.e., snovj electrical characteristics). Under the circumstances the ground 
i truth had to be lim.lted to: 

I a. Density profile of the snow 

I h. Visual inspection of layers in the snow and granularity 

; c. Dielectric constant measurement of a few’ samples, near the top and middle 
I of the snowpaclc 

j d. Wetness measurement of these samples 

I '■ 

Discussion of Experimental Results : The solid line of Fig. 1 shows the backs cat ter ; 

in: dB versus frequency in the range of 5.8 to 8.0 GHz, at the fixed angle of 39® I 
from the vertical, and a snow depth of 70 cm. The electric field vector of the 

I 

horns (20 dB gain each for transmitter and receiver) was parallel to the earth 
surface . 

i The various peaks are produced by layering effects in the snow. The peak 
values are about 25 dB greater than the minima. 

I To then determine what effect frequency deviation, (more specifically, frequency 
laodulation) had on the measurement the signal source v;as reset to the frequency for 
each maximum and minimum, but a modulation of +100 MHz was superimposed for the + 
points; aiTud a modulation of +200 MHz was superimposed for the 0 points. Evidently 
the effect of frequency modulation is merely the averaging of the unmodulated 
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; response over the bc\ndwJdth represented by the riodulat3.on. 

i In Che course of taking ineasui-'ciaeuts , wliethcr at normal incidence or otherx^rLse 

' it was observed Chat the backsca t t er iix g data depended on the observation time and 
i tjemperature . Figures 2> 3 and 4 show the backscatter response at normal incidence, 
i Tlie experimental conditions for the three sets of data v?ere identical except for the 
I time taken and accompanyiiig aixibient temperature. Variation of baclcscatter with snow 

i conditions, at 13.5 GHz and incidence angle of 39" is shovm in Fig. 5. 

i ■ ; 

;■ The decrease in backscatter shortly before noon, was initially observed by us in 

! 1-Iarch 1976 for a snoxvpack 150 cm deep, and has been confirmed at o-ther frequencies 
j by several groups (Currie et al., 1976, Ulaby et al., 1977). Evidently the increase 
I in surface snow wetness by solar heating is the reason for the backscatter decrease, 

! totaling about 15 dB. The air temperature variation' included is useful in 
correlation the phenomenological effects. 

Snow wetness produced during a few hours of sunlight can significantly affect 
the frequency response. March 5, 1977, when the laeasure'uients x^ere taken, Xv'as a 
clear, warm day with a temperature of about 10" C betx/een noon and 16:00, so the 
snow VTas essentially saturated xoritli X’jater in the upper surface. 

Figure 6 shov7s typical, data demonstrating that variation in backscatter as a 
function of frequency can be observed even for other than normal incidence (39") . 
Significant backscatter data has been obtained for incidence angles approaching 45”. 

A continuous variation of incidence angle from the verfical is shoxm 'in Figure 
7 ,: giving the backscatter in relative dB for the fixed frequency of 11.8 GHz (no 
modulation, xv’ith snow depth of 53 cm. The baclcscatter is doxra about 20 dB at the 
angle of 40" from the vertical. 

Theoretical Portion of the Project — K. K. Kei and S. Coen: It is suggested that 

cettain useful properties of the snow-pack may be approximately determined by using 
Inverse scattering theoi^; the inverse problem is that of determining the properties 
of inaccessible snow pack from the back-scattered electric field due to normally 
Incident p.lane waves at varies discrete frequencies. 

Let the back-scattered frequevicy response of the snoxj'-pack be taken as the 
reflection coefficient amplitude jp(o)^)|, i = 1,2,...,N. Neglecting for the moment 
the experimental errors, the question is: how can we learn as much as possible 

about the characteristics of the snow-pack from a set of nunxbers jp(u)^)j, 

1“ 1,2,. ..,N? Alternatively wlrat do the. numbers |p(Wj^) j, 1 = 1,2,...,N tell us 






i about the dielectric constant and conductivity of the snow-pack? The density and 
i wetness of the snow may be then approximately determined by using Weiner’s theory pf 
j dielectric mixtures or other available experimental relations. 

I In what follcras, we show that there is enough information in the numbers 

I |p(Wi)|» 1 = 1,2,...,N, to construct local averages of the snow profile. | 

Hers the experimental errors are neglected and a mathematical model that 
i describes the physical behavior of the snow-pack is chosen. The simplest model 
I that we choose, is that of a plane stratified snow-pack, whose properties are 
I characterized by a linear, isotropic and nonhomogeneous dielectric, medium, as sho\im 
! in Fig. 8a. Maxwell’s equations now show that, for time harmonic dependence 
I the electric field E, satisfies the one dimensional wave equation 

v2 
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{-^-x + k\(x,k)} E(x,k) = 0 


( 1 ) 


3x‘ 
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^diere k = u(ii e )■*■ is the free space wave number and q(x,k) = e(x) - j is the 
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profile function. Here e(x) is the dielectric constant of the snow and cr(x) its 
conductivity. 

Consider two media as shown in Fig. 8b. c. The electric field in the respective 


medium satisfies 


I 


{•^-x + k\ (x,k)} E (x,k) = 0 
^ X n n 


n = 1,2 


( 2 ) 


ind the scattered fields are outgoing. Equation C2) can be reformulated into an 

integral equation form suitable for iterative solutions, 

.1 

■ .h 


6p(k) 6q(x,k) E^(x,k)dx + a(k) 


(3) 


where we have used the relation E 2 - E^^+oE and 6p = P 2 ”Pi were given from 

measurement without error. It may be shown that the term o(h) goes to zero as 
where B ■ 1 is the norm. Tlius, if 6q is sm^ll enough, this term may be 

dropped^ and we obtain a linear approximation 


6p(k) 
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f“ 2 

j 5q(x,k) Ej^(x,k) dx 


( 4 ) 
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viiich is a linear Fredhoira integral equation j, for the 'deterraination of 6q from the 
given functions 5p and E^. Once 5q is found, the profile is given by 

“ %+‘5q. Ucrwever, equation (4) is not exact, so the profile will not in 
general, satisfy the equation (2). To overcome this, v?e solve the equation by 
iterations, until convergence occurs. In practice fev^ iterations will suffice. 

1 la practice, the amplitude of the reflection coefficient is much to be 

preferred than its phase, and equation (4) can be rewritten into the form 
.h ■ 
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^ 6e(x) F^(x)dx 


1 * 1.2, 


.W 


(5) 
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where = log 


PjCki) 


Pl(lc^) 


5e 


and F^(x) 


Re{-r^ — TT^ E, (k. ,x)}; here 
2jp,(k.) 1 X 


the frequency k is discretized and the medi-um is assuaied lossless. The approach to 
be discussed is however applicable to lossy medium as well. 

We now shcA .7 how to solve Eq. (5) by the Backus and Gilbert approach. Mult5.ply 
the equation by (5) and sum over i; this results in 


( 6c ), 
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5e(x) A(C.X)dx 


where 


N 

V* 


<«e> = ^.a) g. 

^ • i=l 


(6) 


(7) 


Is the local, average of 6 g» and 


N 


A(^,x) - aAo f .(^) 

i=l 


( 8 ) 


|lsi the averaging kernel. 

The coefficients 0-^iO • 1 “ 1,2,. ..,N are determined, such that the averaging 
kernel resembles the Dirac-Delta function most closely. One way of doing this is by 
minimizing the quadratic form 


i: 


(C-x)^A^(C.x)dx 


(9) 
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subject to . 

A(C,x)dx = 1 (10) 

I for each C» this minimization problem may be reduced to the linear equation upon 
I introducing a Lagrange multiplier. 

I Once the coefficients are determined, the local average (5e). is 

i ^ 

I computed via equation (7) , 

In Figs. 9 and 10 we px'esent results for 3 different snow profiles; the results 
I shown in Figs 9a and 9b have been obtained with data at 3 frequencies only, 

whereas, the results sho\'ra in Fig. 10 have been obtained with data at 7 frequencies. 
The approximate results obtained from our analysts are in good agreement with the 
exact results for smooth profiles, as is evident from Figs. 9a and 9b. The poor 
resolution near the edges of the ice layer (Fig. 10) may be improved by using more 
data. It will be noted that the computation cost is about 5 dollars per profile 
on CDC 7600. 

We are currently investigating the influence of e:-:perimental errors on the 
resulting snow profiles. 

Wc find that there is a tradeoff between resolution and uncertainty, which is 
Ini agreement with the Backus and Gilbert theory. 
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